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ABSTRACT 

We measured stellar velocity dispersions a and derived dynamical masses of 9 massive {M « 10^"'^ 
Mq) early-type galaxies (ETG) from the GMASS sample at redshift 1.4 < z < 2.0. The a are 
based on individual spectra for two galaxies at 2 w 1.4 and on a stacked spectrum for 7 galaxies 
with 1.6 < z < 2.0, with 202-h of exposure at the ESO Very Large Telescope. We constructed 
detailed axisymmetric dynamical models for the objects, based on the Jeans equations, taking the 
observed surface brightness (from deep HST/ACS observations), PSF and slit effects into account. 
Our dynamical masses Mjcans agree within < 30% with virial estimates M^ir = 5 x R^a^ /G, although 
the latter tend to be smaller. Our Mjeans also agrees within a factor < 2 with the Mpop previously 
derived using stellar population models and 11 bands photometry. This confirms that the galaxies are 
intrinsically massive. The inferred mass-to- light ratios {M/L)u in the very age-sensitive rest frame 
[/-band are consistent with passive evolution in the past ^ 1 Gyr (formation redshift 2;^ ~ 3). A 
'bottom-light' stellar Initial Mass Function (IMF) appears to be required to ensure close agreement 
between Mjeans and Mpop at z ^ 2, as it does at z ~ 0. The GMASS ETGs are on average more 
dense than their local counterpart. However a few percent of local ETGs of similar dynamical masses 
also have comparable a and mass surface density E50 inside Re- 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: evolution — galaxies: formation 
— galaxies: high-redshift 



1. INTRODUCTION 

In the hierarchical galaxy formation paradigm, where 
galaxies are assembled by the merging of multiple build- 
ing blocks in a universe dominated by dark matter (e.g. 
Springel et al. 2005), the most massive early-type galax- 
ies (ETGs) are assembled last. However observations in 
the local universe and at high redshift seem to converge 
towards a 'downsizing' mechanism for ETGs formation in 
which the the stars of the most massive systems formed 
at the highest redshifts (z > 3), while the stars in the 
smaller ones were produced over more extended periods 
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of time (Cowie et al. 1996; Heavens et al. 2004; Thomas 
et al. 2005; Treu et al. 2005; Renzini 2006). A way to 
reconcile the apparent contradiction between these two 
pictures is to assume that the stars in the massive sys- 
tems formed via efficient star formation processes at high 
redshift and were later assembled into larger systems via 
mostly collisionless mergers (De Lucia et al. 2006; Khoch- 
far & Silk 2006; Naab et al. 2009). 

An important test for this scenarios is constituted by 
the mass and size distribution of ETG at z > 1. Contrary 
to the expectations the most massive ones appear to be 
already in place (Cimatti et al. 2004, 2006; Glazebrook 
et al. 2004; Scarlata et al. 2007) but have much smaller 
sizes than their local counterparts (Daddi et al. 2005; di 
Serego AHghieri et al. 2005; Trujillo et al. 2006). This 
suggests that they might not be the direct precursors of 
present-day ETGs and mergers must play a role in their 
evolution (Trujillo et al. 2007, 2009; Zirm et al. 2007; 
Longhetti et al. 2007; Toft et al. 2007; Cimatti et al. 
2008; van Dokkum et al. 2008; van der Wei et al. 2008; 
Buitrago et al. 2008; Bernardi 2009). 

There are however concerns affecting the compactness 
determinations due to possible observational biases af- 
fecting either the mass or size estimate (due to the sur- 
face brightness dimming, the presence of AGNs or nu- 
clear starbursts) of galaxies at high redshift. Here we 
try to to address these concerns by measuring the veloc- 
ity dispersion a of the stars, related to the density, and 
deriving masses via dynamical models of mass-selected 
ETGs at 1.4 < z < 2.0. We assumed a flat Universe 
with Ha = 70 km s^^ Mpc"\ f7,„ ^ 0.3, f^A = 0.7. 

2. VELOCITY DISPERSION DETERMINATION 
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2.1. Spectroscopic and photometric data 

The sample under exam comes from the Galaxy Mass 
Assembly ultra-deep Spectroscopic Survey (GMASS)^'^ 
within the redshift range 1.4 ^ z < 2.0 (Cimatti et al. 
2008, hereafter COS). It was flux-selected at 4.5 /xm using 
the Great Observatories Origin Deep Survey GOODS- 
South public image taken with IRAC on the Spitzer 
Space Telescope (Dickinson et al., in preparation). 

The GMASS optical multi-slit spectroscopy used here 
was obtained with the ESO VLT + F0RS2 (MXU mode) 
in the wavelength range 600-1000 nm with the grism 
3001, using very long integration times of up to 32 h 
per spectroscopic mask, and with a slit width of 1 arc- 
sec. We adopted as instrumental resolution the mean 
Cinstr = 130 ± 21 km s~^ of the values derived from sky 
emission lines and a star, where the error is half the 
difference between the two determinations. This con- 
servative error also accounts for the small dependence 
of the resolution with wavelength. We also use public 
HST/ACS/F850LP (z-band) photometry from GOODS- 
South (Giavalisco et al. 2004). 

2.2. Library of stellar templates 

The F0RS2 observations span a rest-frame UV wave- 
length range of 230-385 nm at the mean redshift z ^ 1.6 
of the GMASS sample. To measure stellar kinematics 
we need stellar templates in the UV and we cannot use 
the extensive ground-based stellar libraries. Moreover 
no empirical UV library span the full required spectral 
range. 

For this we use synthetic libraries, which now can re- 
produce spectra of real stars remarkably well (Munari 
et al. 2005; Martins & Coelho 2007). The mismatch 
in minor spectral features is not critical when working 
with low-5'/A^ spectra dominated by systematics. Here 
we selected as templates a subset of 33 models from 
the high-resolution R = 20, 000 synthetic spectral li- 
brary-'^'* by Munari et al. (2005) spanning a wide range of 
temperatures 3500 < T < 10,000 and surface gravities 
< log 5 < 5, at solar metallicity and abundance. 

2.3. Individual spectra at z lA 

Some sharp absorptions are required for a reliable kine- 
matics extraction. In the restframe wavelength range of 
interest (250-400 nm) the only significant ones are the 
Mg II doublet (280 nm), Mg I (285 nm) on the blue side, 
the Ca II H and K (~ 395 nm) absorptions and a blend 
of Fe I and Mg I (384 nm) on the red side. 

At z > 1.5 the red spectral features fall outside our 
observed red range of 1000 nm and the kinematics re- 
lies on the Mg II doublet and Mg I absorptions. Suffi- 
cient S/N is required for robust measurements at these 
redshifts. For this we could only measure reliable a for 
two z w 1.4 individual galaxies, where the red features 
could be included in the fit. For all a determinations 
we used the Penalized Pixel-Fitting mcthod^^ (pPXF; 
Cappellari & Emsellem 2004) with the 33 templates of 
Munari et al. (2005), including additive polynomials, to 
correct for residual template mismatch or sky subtrac- 
tion errors, and multiplicative polynomials, to correct 

http:/ /www. arcetri.astro.it/~cimatti/gmass/gmass. html 
1* http://archives.pd.astro.it/2500-10500/ 

http:/ /www-astro. physics. ox. ac.uk/~mxc/idl/ 
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Fig. 1.— Kinematics extraction for GMASS 2470 (z Ri 1.4). 
The panels show the pPXF fits to the full spectral range (Top 
Panel), the blue (Middle Panel) and the red one (Bottom Panel) 
respectively. In each panel the black line is the observed spectrum, 
the red one is the best fitting template and the green diamonds are 
the residuals (arbitrarily shifted). The two vertical dashed lines 
indicate the spectral region excluded from the fit due to high noise 
due to sky lines. 

possible spectral calibration errors. We verified the sta- 
bility of our results with different degrees between 1-4 
for the two sets of polynomials. In the fits 4-8 of the 
33 templates were selected by pPXF to reproduce the 
spectrum. 

The spectrum with the highest mean S/N « 8 is 
GMASS 2470 (table 1 of C08). We measured CTobs (before 
correcting for instrumental resolution) for three wave- 
length ranges (Fig. 1): (i) The full range (255-405 nm): 
(Jobs = 192 ± 13; (ii) The blue range (264-297 nm): 
(Jobs = 198 ± 28; (iii) The red range (371-408 nm): 
Cobs = 159 ± 18. The three results are consistent within 
the relative error bars, giving confidence in the adopted 
approach. We found in general no trend with wavelength 
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TABLE 1 

Sample of GMASS passive early type galaxies and measured parameters 



ID 

(1) 


z 

(2) 


AV 
(km s~ 
(3) 


•^prcd 

^) (km s~ 
(4) 


1) (km s-l 
(5) 


Act* 
) (km s-1) 
(6) 


(7) 


He 
(arcsec 
(8) 


logLu (Ai/L) Jeans 
:) (kpc) (Lq^) (Mq/Lq^) 

(9) (10) (11) 


log Mjeans 

(12) 


log Mvir 

(Mo) 
(13) 


log Mpop 
(Mq) 
(14) 


0472 


1.9077 


36 


191 






2.7 


O.06 


0.54 


11.04 


0.38±0.09* 


10.64" 


10.53'' 


10.49 


0996 


1.3844 


30 


98 


— 


— 


2.4 


0.13 


1.06 


10.36 




— 


— 


10.16 


1498 


1.8491 


16 


157 


— 


— 


2.5 


0.14 


1.18 


10.94 


0.98±0.22* 


10.93* 


10.83" 


10.61 


2111 


1.6102 


19 


185 


— 


— 


4.0 


0.09 


0.80 


10.85 


0.79±0.18^ 


10.75'' 


10.68" 


10.61 


2148 


1.6118 


23 


248 






6.5 


0.14 


1.22 


11.00 


0.89±0.20^ 


10.95'^ 


10.84" 


11.02 


2196 


1.6063 


28 


180 






3.1 


0.17 


1.40 


10.89 


0.92±0.21^ 


10. 85=^ 


10.89" 


10.79 


2239 


1.4149 


15 


113 


Ill 


35 


4.5 


0.25 


2.09 


10.60 


1.23±0.78 


10.69 


10.52 


10.54 


2286 


1.6020 


29 


135 






3.2 


0.18 


1.48 


10.72 


1.73±0.39^ 


10.95=' 


10.91" 


10.56 


2355 


1.6097 


32 


127 






2.2 


0.12 


1.00 


10.78 








10.36 


2361 


1.6096 


18 


197 






4.1 


0.15 


1.26 


10.86 


0.93±0.21^ 


10.83" 


10.85" 


10.83 


2470 


1.4149 


11 


157 


141 


26 


7.6 


0.18 


1.53 


10.92 


0.66±0.24 


10.74 


10.61 


10.71 


2543 


1.6149 


54 


141 






1.8 


0.22 


1.88 


10.60 








10.69 


2559 


1.9816 


28 


147 






2.4 


0.19 


1.61 


10.91 








10.67 


New Stack l.( 


3 < 2 < 2.0 


15 


205'' 


202 


23 


8.0 




1.16*^ 


10.93" 


0.93'= 


10.88" 


10.82" 


10.85" 


Public Stack 


all 


13 


175^= 


< 214 




8.7 




1.37^= 


10.88<: 








10.76== 



Note. — Column (1): GMASS ID from C08. Column (2): Redshift measured with pPXF. Column (3); \a error in the velocity 
alignment. Column (4): Virial prediction for the velocity dispersion a% = GMpop/{5Re), corrected to a 1x1 arcsec^ aperture. Column 

(5): Measured galaxy velocity dispersion. Column (6): Error on cr*. Column (7): S/N per 60 km spectral pixel, computed from the 
pPXF fit residuals to the GMASS spectrum. Column (8): Circularized Rc from the PSF-deconvolved MCE model. Column (9): Rc in 
kpc. Column (10): Total ?7-band luminosity from the MGE model. Column (11): C/-band mass-to-light ratio from the Jeans dynamical 
model. Column (12): Total mass from the Jeans model. Column (13): Virial estimate of the total mass. Column (14): Stellar population 
estimate of the mass from C08, using the models of Maraston (2005) normalized for a Chabrier (2003) IMF. 

"The spectrum of this galaxy was included in the New Stack. 
These values were computed by adopting for the galaxy the cr* of 
the stacked spectrum of Section 2.5. 

^Weighted mean {u) = Ej [«j (^/Af)!]/ Ej {S/N)^j of the quanti- 
ties Uj for the 7 galaxies included in the New Stack (Section 2.5). 

"^Weighted mean of the quantities for all 13 galaxies in the Public 
Stack (Section 2.4). 



SO wc adopt as standard value the one measured for the 
full speetral range, which as expected has smaller errors. 
The galaxy stellar dispersion ct* is: 



(1) 



introduce an additional Gaussian velocity broadening 

attack = AV^Az c/(l + z) (3) 

in the coadded spectrum. Assmning all broadening func- 
tions to be Gaussian, the dispersion of the individual 
galaxies could be recovered using 



The values and errors for this galaxy and for GMASS 
2239 are given in Table 1. 

2.4. Public GMASS Stacked spectrum 

We applied the same approach of Section 2.3 to mea- 
sure (T* from the public GMASS spectrum, obtained by 
coadding the individual normalized spectra of 13 ETG 
within 1.4 < 2; < 2.0, for an unprecedented total integra- 
tion time of 480-h and a mean S/N « 9 (see fig. 4 of COS). 
We derived (Tobs = 287 ± 20 km s^^. Our value and er- 
ror agree with the determination (Tobs = a/ST^ -f 267^ = 
273 ± 20 km (with high-resolution UV stars) per- 
formed on the same spectrum by Cenarro & TVujillo 
(2009). Those authors used empirical stellar template 
spectra, so this agreement validates our approach of us- 
ing synthetic templates in the UV. 

An additional step is needed to estimate the character- 
istic a* of the galaxies in the stack. In fact the spectra 
had to be shifted to the rest frame wavelength before 
coaddition. The measured redshift can be written as 



obs 



stack ■ 



(4) 



1 + Z = (1 + Ztruc) X (l + V/c), 



(2) 



where ztmc is the true galaxy redshift, V is the velocity 
shift due to an error in z and c is the speed of light. 
If the galaxies had all identical spectra and the redshift 
errors Az were normally distributed, the stacking would 



In the case of the public GMASS spectrum, which was 
not intended for kinematics measurements, the individ- 
ual z were measured via cross correlation and the small- 
est (Tstack ~ 141 km s~^ (M. Mignoli private communi- 
cation). One can then only derive an upper limit to the 
typical dispersion in the stack tr* < 214 km s^^. This 
limit is smaller than the cr* derived by Cenarro & Trujillo 
(2009) as they incorrectly assumed (Tstack to be negligible. 

2.5. New stacked spectrum oi 1.6 < z < 2.0 

We re-measured z of ah 13 GMASS galaxies with pPXF 
and give redshifts and errors in Table 1. The new aver- 
age velocity-error becomes CTgtack ~ 30 km s~^, which is 
negligible with respect to the expected dispersions. We 
verified the reliability of our errors by measuring z of the 
individual exposures of the same galaxy. 

After excluding the two spectra of Section 2.3, to max- 
imize the S/N we constructed a stacked spectrum from 
the 7 remaining GMASS spectra with S/N > 2.5. We 
normalized the spectra in the 260-310 nm wavelength 
range before coaddition. not to bias the kinematics to- 
wards the brightest galaxies. The measured (T^ = 202±23 
km s~^ for the stack, corrected with equation (4), agrees 
with the weighted average (fTpred) = 205 km s~^ of the 
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stack of 7 galoxles 1.6 < z < 2.0 (202-h) 




Fig. 2. — Kinematics from stacked spectrum of Section 2.5. In 
each panel the black line is the observed spectrum, the red one is 
the best fitting template and the green diamonds are the residuals 
(arbitrarily shifted). The blue crosses indicate pixels automatically 
excluded from the fit. The solid blue line indicates the estimated 
1(7 noise. 
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Fig. 3. — MGE models for ten GMASS galaxies for which we 
constructed dynamical models. The contours of the observed 
HST/ACS/F850LP surface brightness are overlaid to an MGE 
model of their surface brightness, convolved with the ACS PSF. 
Contours are spaced in 0.5 mag arcsec"^ intervals. 

virial predictions (Tabic 1) for the galaxies in the stack. 
Although we do not trust the individual measured tr^ val- 
ues for each \ow-S/N spectrum in the stack, and do not 
give them in this paper, they are also not inconsistent 
with tJpicd and span the same range of values. The cr^ 
of the two individual galaxies of Section 2.3 also agrees 
with CTprcd- In Table 1 and in what follows we adopt the 
CT* from the stack as representative of the cr* of each of 
the 7 galaxies in the stack. This is not correct for each 
individual case, but only in an average sense. 

3. DYNAMICAL MODELS 
3.1. Jeans modeling 



The (T+ we measured for the GMASS galaxies in Sec- 
tion 2 can be used to determine their dynamical masses. 
As the galaxies have half-light radii i?e ^ 0'.'25, while the 
spectra aperture and seeing have size of ~ 1" one may 
need significant corrections to the virial formalism used 
at low redshift to measure masses. One may estimate cor- 
rections using spherical Sersic dynamical models based 
on the Jeans equations (di Serego Alighieri et al. 2005). 
However these models cannot describe well all ETGs, 
especially when they have disks and may rotate signif- 
icantly. For this reason van der Marel & van Dokkum 
(2007) and van der Wei & van der Marel (2008) used ax- 
isymmetric Jeans dynamical models of individual galax- 
ies to take the surface brightness and possible rotation, 
as well as PSF and aperture, directly into account when 
measuring masses at high redshift. This is the approach 
we also use here. 

We adopt a Multi-Gaussian Expansion (MGE) 
(Emsellem et al. 1994) to parameterize the 
HST/ACS/F850LP (z-band) surface brightness of 
the GMASS galaxies (Fig. 3), while taking the ACS 
PSF into account, using the software^^ of Cappellari 
(2002). The following expression was used to if -correct 
the MGE parameters from observed countrate {Cz), 
in counts s"""^ per ACS pixel, into a restframe Johnson 
?7-band surface brightness in mag arcsec"^ 



M(7 = -2.5 log 



C'z X /850 X (1 + z)' 



(5) 



fu{AOV) X p2 

Here fs5o = 1-51 x 10^^^ ergs s"^ cm"^ is the lat- 
est inverse sensitivity of the F850LP filter,^^ fuiAOV) = 
4.28 X 10^^ ergs s^^ cm^^ A^^ is the zero point of the 
Johnson [/-band, p — O'.'OS is the dithered pixels size 
of the GOODS images." We include both the (1 + z)^ 
bolometric dimming of the surface brightness and a fac- 
tor (1 + z) due to the redshifting of the bandwidth. The 
formula is accurate at z w 1.4, where the ACS/F850LP 
band is de-redshifted into the U -band. At larger redshifts 
we applied a small extra JsT-correction inferred from the 
stacked GMASS spectrum. 

For each redshift we placed the models at the corre- 
sponding angular diameter distance Da- We computed 
a prediction for the velocity second moment (Vj.^g = 
+ cr^) inside a 1" square aperture, with a 1" seeing 
FWHM, assuming semi-isotropy (/3^ ~ 0) and axisym- 
metry, for a constant {M/ L)u — 1, using equation (28) 
of the Jeans Anisotropic MGE (JAM)^^ method of Cap- 
pellari (2008). We assumed an intermediate inclination 
i = 60° for all galaxies, but the results do not change 
more than 5% for an edge on inclination (i = 90°). 
The dynamical M/L of each galaxy is then given by 
Jeans — (f^ * /K-ms)^ (Table 1). The (M/L) 

Jeans 

decreases by 5% by assuming in the models the largest 
radial anisotropy (3z — 0.5 observed in nearby galaxies. 

Not all 7 galaxies included in the stack are expected to 
have the same tr* — 202 km s~^ we measured. Some can 
be higher and some lower than this average. If the virial 
predictions CTprod were correct, the fact that generally 
fprcd < 202 km s^^ suggests the quoted masses and M/L 
are mostly overestimated. 

http:/ /www. stsci.edu/hst/acs/analysis/zeropoints 

http;/ /archive. stsci.edu/pub/hlsp/goods/v2/h_goods_v2. 0_rdm.html 
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3.2. Virial and population masses 

In the top panel of Fig. 4 we compare the dynam- 
ical mass obtained from the JAM models Afjoans = 
L[/X (M/L) Jeans to the virial massMvir = 5.0x i?cO'c^/G, 
where the scaling factor was calibrated using dynami- 
cal models and integral-field data of local ETGs, for CTo 
measured within IRc (Cappellari et al. 2006). We es- 
timate (Te by increasing cr^ from the measured 1" x 1" 
aperture to a IRc circular aperture using equation (1) 
of Cappellari et al. (2006). Our values were deter- 
mined in a non-parametric way from circularized MGE 

models {{crj,qj} ^ Wi^j^jA}) which preserve the lu- 
minosity and peak surface brightness of each Gaussian. 
With constant ellipticity this corresponds to the circu- 
larized radius i?o — Vab of the ellipse enclosing half of 
the analytically-derived MGE galaxy light. Our values 
agree (except for GMASS 2196) with the determination 
via Sersic profiles fits of COS within their quoted 20% 
errors (estimated via simulations). 

There is a general agreement between the JAM and 
virial estimate, but in median the latter is ^ 30% lower. 
As both values are based on the same ct*, the difference 
must be attributed to an underestimation of i?e and/or 
to non-homology in the profiles. This may be due to the 
low S/N caused by cosmological surface brightness dim- 
ming (Mancini et al. 2009). The JAM approach has the 
important advantage over the virial one that it robustly 
recovers the M/L even when non-homology is important 
or the outer parts of the profiles are lost in the noise. 
Considering a test model with an I{R) cx exp(— fci?^/"*) 
surface brightness profile truncated at li?o, we still re- 
covered the true M/L to 1% with JAM, but the M/L 
was underestimated by 26% with the virial approach. 

In the second panel of Fig. 4 we compare Mjeans to 
the mass determination Mpop based on stellar popula- 
tion models and 11 photometric bands of COS. The val- 
ues are in agreement within the rather large uncertainty 
in both quantities. The agreement may improve when 
considering the possible underestimation of cr^, for some 
galaxies in the stack. This shows that mass errors are 
< 2x when detailed photometric information is avail- 
able. It also confirms the result of COS that ETGs at 
z ^ 2 are consistent with a passive evolution in the past 
t ^ 1 Gyr and indicates a formation redshift z ~ 3. 
Any significant star formation activity would have dra- 
matically lowered the dynamical {M/L)u which scales 
linearly with time in the age-sensitive U-hand. This is 
in agreement and extends to z ~ 2 previous dynamical 
studies oi AI/L evolution based on the the Fundamental 
Plane at z ~ 1 (van Dokkum & Stanford 2003; Gebhardt 
et al. 2003; van de Ven et al. 2003; van der Wei et al. 
2004; Treu et al. 2005; di Serego Alighieri et al. 2005, 
2006; J0rgensen et al. 2006). 

The Mpop values are based on the Chabrier (2003) 
Initial Mass Function (IMF). Adopting a Salpeter IMF 
would increase Mpop by 70%, making in most cases 
Mjeans < Mpop for thcsc high redshift galaxies. Similarly 
with a straight Salpeter IMF the (M/L)pop ratio of local 
ETGs would be about twice the value derived from dy- 
namical modeling (Renzini 2005; Cappellari et al. 2006), 
which requires instead a bottom-light IMF such as in the 
case of Kroupa (2001) or Chabrier's IMFs. Therefore, it 
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M.ecns (Mo) 

Fig. 4. — Top Panel: Comparison between mass determinations 
via dynamical models Mjeans and virial masses Mvir- Blue dia- 
monds are galaxies with individually measured tr^ (Section 2.3), 
while red circles are galaxies for which we assumed the (t^. of the 
stacked spectrum (Section 2.5). The dashed line indicates equality, 
while the dotted lines correspond to a factor 2x difference. Second 
Panel: As in the top panel, for a comparison with stellar popu- 
lation masses Mpop. The error bars in the latter span the ranges 
of estimates using the three different population codes presented 
in COS. The symbols correspond to the Maraston (2005) models. 
Third Panel: Comparison between the GMASS (Tc, the values for 
the Coma sample (black filled circles) and dynamical models of lo- 
cal ETGs (black open squares). The large red open circle indicates 
the weighted average value for the seven galaxies in the New Stack 
(Table 1). The dash-dotted lines are virial predictions of mass sur- 
face density S50 inside Re- Bottom Panel: As in the third panel, 
for E50. We adopted errors of 30% on Re- The dash-dotted lines 
are virial predictions of (Tg. The solid green lines are the values 
and errors from Shen et al. (2003). The magenta arrow represents 
a 30% decrease of cr* . 
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appears that the dynamical modeling of both low redshift 
and high redshift ETGs requires a bottom-light IMF. 

4. DISCUSSION 

We have mcasiired the stellar velocity dispersion a, 
from individual and stacked spectra, and have con- 
structed detailed dynamical models, of 9 early type 
galaxies (ETGs) from the GMASS sample (C08)' in the 
redshift range 1.4 < 2; ^ 2.0. The agreement between 
the dynamical masses and the ones previously derived 
via population models by COS indicates that an overes- 
tiniation of the mass can not explain the high density 
discovered by previous works. 

If high-z ETG are indeed denser that local ones, they 
should have a higher a and surface mass density S50 = 
Mjeans/ (27ri?c^) withiu i?c at given dynamical mass (Toft 
et al. 2007; van Dokkum et al. 2008). To test this fact, 
in the bottom two panels of Fig. 4, we compare the mea- 
surements for our GMASS galaxies to a sample of ETGs 
in the Coma cluster (J0rgensen 1999; J0rgensen et al. 
2006) and to dynamical models of local ETGs (Cappel- 
lari et al. 2006) , which use the same modeling technique 
as this paper. Wc also compare with the density derived 
on SDSS galaxies by Shen ct al. (2003), increased by 30% 
to account for the fact that the population masses using 
a Kroupa IMF on average underestimate the dynamical 
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mass of massive early-type galaxies (e.g. fig. 17 of Cap- 
pellari et al. 2006). We find that our two z w 1.4 galaxies 
have a and S50 consistent with the ones of local ETG (as 
shown in COB). However the galaxies in the stacked spec- 
trum at 1.6 < 2: < 2.0 have on average the a and E50 of 
the most dense local ETGs. 

This paper illustrates the limits of what can be 
achieved on the study of the dynamics of ETG with the 
current generation of telescopes. It emphasize the use- 
fulness of stacking technique to infer the dynamics of 
selected classes of galaxies. Much progress along these 
lines could be obtained with massively multi-object spec- 
trographs on the future generations of 30 40 m tele- 
scopes like the E-ELT. Access to an atmosphere-free 
near-infrared wavelength range, as soon available on 
James Webb Space Telescope, would dramatically im- 
prove the kinematics determination in ETGs at z > 2 by 
bringing the rich set of optical absorption lines into the 
observable domain. 
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